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Abstract 


In this note we study an ideal transmission line, which is described 
by the telegraph’s equation with variable coefficients, from the view- 
points of numerical analysis and control theory. Because the semi- 
discretized systems of original system are not suitable to discusses the 
uniform exponential stability, we perform a similar transform to the 
continuous system and obtain an intermediate system, which is con- 
veniently studied. Firstly, we investigate uniform exponential stability 
for the intermediate system by a so called average central-difference 
semi-discretization scheme and direct Lyapunov function method. The 
proof is parallel to that of the continuous case. Secondly, the stabil- 
ity and consistence of numerical approximating scheme are presented 
by the Trotter-Kato Theorem. At last, a semi-discretization scheme of 
the original system is proposed by the inverse transform. All results 
for the intermediate system are then translated into the original sys- 
tem. As the byproduct of the main results, the uniform observabil- 
ity problem is also addressed. Furthermore, the effectiveness of the 
numerical approximating algorithms is verified by numerical simulations. 
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1 Introduction 


Space semi-discretization, which transforms PDEs into ODEs, is a first natu- 
ral step in the process of numerical discretization of PDEs. Many researchers 
addressed this problem from different viewpoints in the past decades. Math- 
ematicians mainly focused on stability and consistency of the approximating 
algorithm and designed many profound discrete schemes [1]-[3]. Physicists 
aimed at preserving structural invariants and structural properties of the 
continuous system. Several recent works tackled this challenge for infinite- 
dimensional port-Hamiltonian systems. The Dirac structure of the system is 
preserved by mixed finite element or finite difference methods [4]-[6]. 

However, when preserving control properties, such as passivity, exponential 
stability and observability etc, are considered, detailed and elegant analysis 
is further needed on the basis of stability and consistency of the approximat- 
ing algorithm [7]-[11]. Banks, Ito and Wang firstly pointed out that the wave 
equation with boundary damping didn’t inherit the exponential stability of 
the continuous one when the wave equation is discretized by the classical finite 
difference and finite element schemes in [12]. In the same time, Glowinski, 
Li and Lions showed that the exact controllability property wasn’t preserved 
for some discretization process in [13]. Lately, Infante and Zuazua discussed 
the problem of boundary observability of wave equation, i.e., the problem of 
whether the total energy of solutions can be estimated uniformly in terms of 
the energy concentrated on the boundary as the net-spacing approaching zero. 
They derived negative answer due to the existence of high frequency spuri- 
ous solutions for both finite-difference and finite-element semi-discretizations. 
A uniform bound was obtained in a subspace of solutions generated by the 
low frequencies of the discrete system [14]. This method is in fact a filtering 
strategy which has been widely used in control theory [15]-[16]. Zuazua gave 
comprehensive survey on the observation and control of waves approximated 
by finite difference methods in [17]. 

To overcome the obstacle caused by high frequency spurious modes, several 
solutions were proposed. One suggestion is to use mixed finite element method 
[4]-[6] to obtain the uniform controllability of the conserved wave equation [18]. 
Some other remedies of damping out high frequencies include Tychonoff reg- 
ularization [13], two-grid algorithms [19] and non-uniform numerical meshes 
[20]. Introduction of a vanishing viscosity term at whole domain of the spatial 
variables is another popular approach [21]-[23]. Very recently, Liu and Guo 
introduced an average operator for the time derivative of classical finite dif- 
ference for the wave equation with boundary damping and showed that the 
scheme uniformly preserves the exponential decay of the continuous system 
[25]. Using the same idea, Xu, Guo and Zheng proposed two completely differ- 
ent finite-difference schemes for uniform exponential approximations of wave 
equation with local viscosity damping [26]-[28]. Guo and Zheng et al. gener- 
alized the results of [25]-[28] to the coupled heat-wave system [29] and the 
Schrödinger equation on L?(0,1) space [30]. 
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Inspired by above works, we shall study the transmission line with vary- 
ing capacity and inductance, which is described by the telegraph’s equation 
with variable coefficients, from the viewpoints of numerical analysis and con- 
trol theory in this note. It is well known that in engineering applications, there 
are three common physical systems that carry waves: the electrical transmis- 
sion line, the flexible string, and the compressible fluid [4], [31], [32, Section 
7.1]. Since the electrical transmission line is essentially the wave, the effective 
methods of [25]-[29] can be employed to study it numerically. But the vari- 
able coefficients appearing in the space derivatives of PDEs bring some new 
troubles in construing suitable numerical approximating algorithm to preserve 
exponential stability. To the best of our knowledge, classical finite-difference 
method for the spatial semi-discretization of the telegraph’s equation of this 
note hardly preserves the passivity, not even the exponential stability of the 
continuous system. 

To overcome these difficulties we first perform a similar transform to the 
original system and obtain an intermediate system, which is easily studied. 
A semi-discretization scheme, which is called average central-difference (see 
Section 4), is proposed for the intermediate system. We then investigate the 
uniform exponential stability of the discrete systems utilizing the method par- 
alleling to that of the continuous case. We also reveal how the mechanism 
of boundary feedback control uniformly exponentially stabilizes the discrete 
system (see Remark 2). Secondly, the stability and consistence of numerical 
approximating algorithm are presented by the Trotter-Kato Theorem. At last, 
a discretization scheme of the original system is proposed by another similar 
transform. All results on the intermediate system are then translated into the 
original system. As the byproduct of the main results, the uniform observabil- 
ity is also presented. Furthermore, to show the effectiveness of the numerical 
approximating algorithms, we perform several numerical simulations. 

Thus, the contribution of this work is four-fold: 


e Give a new method to study numerical solution of PDEs with variable 
coefficients. In the meanwhile, explain how the feedback control uniformly 
exponentially stabilize the discrete system. 

e Extend the results of uniform exponential stability of [21] and [25]-[28], and 
uniform observability of [14] and [17] from simple models to more complex 
system. 

e Provide demonstration to study uniform exponential stability of other com- 
plex models described by PDEs such as wave-wave coupled equations, beam 
equations and so on. 

e The results of uniform exponential stability and uniform observability have 
potential applications in uniform controllability, the approximation of the 
control problem and state reconstruction etc. 


The structure of this paper is as follows. In section 2, the transmission 
line system is introduced and some results about exponential stability and 
exact observability of the intermediate system are presented. In section 3, the 
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uniform exponential stability and uniform observability of discrete systems are 
obtained. In section 4, the stability and consistence of numerical approximating 
algorithm is derived by the Trotter-Kato Theorem. In section 5, a discretization 
scheme of the original system is proposed by similar transform. The uniform 
exponential stability of the original system is proved and convergence analysis 
of the discretization scheme is made based on the result of section 4. In section 
6, numerical simulations are performed. In section 7, we give some concluding 
remarks. 


2 Results of continuous model 


In this section, we introduce the model discussed in this paper and present 
some known results about the exponential stability of continuous systems. 
Consider the transmission line on the spatial interval [0, 1]: 


£02) =-2 (#2), ce (0,0, 
t 


ð 
HOt, 2) = -i (43) ,t>0, (1) 
V(t,0)=0, V(@,1)= RI(t,1), R>0 


1 
Q(0, £) = Q(x), (0, x) = ¢° (x). 


Here Q(t, x) is the charge at position x € [0,1] and time t > 0, and ¢(t, x) is 
the magnetic flux at position x and time t. C(x) is the distributed capacity 
and L(x) is the distributed inductance. The voltage and the current are given 
by V = Q/C and I = @/L, respectively. V(t,1) = GI(t,1) is the boundary 
feedback and R is the feedback gain constant. (Q°(-), 6°(-)) € [L£7(0,1)]? are 
initial configuration of the transmission line model. The energy of this system 


is given by i 
1 f lotz)? | Q, x)|? 
E(t) = 5 [ ie + C(x) da. (2) 
It follows from the exercises 7.1 and 9.1 of [32] that the system (1) is expo- 
nentially stable with respect to the energy E(t). In order to study the uniform 
exponential stability in a convenient manner, we begin from the following 
system on V and T: 


C(x)V;(t, x) = —I, (t,x), x € (0,1), 

L(x)I,(t,2) = —V,z (t,x), t > 0, (3) 
0 RI(t,1), 

Oa) = I(x). 


aE aa 
II 


The energy of system (3) is also given by 


E(t) = 5 | C(x)|V(t,2)/? + L(a)|L(t, x) Pde (4) 
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and has different expression in light of V = Q/C and I = ¢/L. To give some 
clues in discrete case, we apply the method of Lyapunov function to obtain the 
exponential stability of the continuous system (3). For this purpose, we assume 
that the capacity function C(x) and the inductance function L(x) satisfy: 

Hı: C(x) > 0, L(x) > 0, Vx € [0,1]; 

Hə: C(x), L(x) € C*[0,1] and C(x) < K, L(x) < K for some positive 
constant K; 

H: C'(x) > 0, L'(x) > 0, Vx € [0,1]. 

It should point out that H3 is not applied in the proof of the main result, 
see for instance Theorem 3. It is only used in the method of Lyapunov function 
to verify the exponential stability of the continuous system. If one applies the 
method of [32, Lemma 9.1.3], H3 is also useless. 


Theorem 1 Under the conditions H,-H3, for any V°(-), [°(-) € L?(0, 1], there exist 
two constants M and w such that the energy of the solution to the system (3) satisfies 


E(t) < Me “'E(0). (5) 


Proof: It is easy to see that 


< Blt) = 1 veantornta ei 
= —|V(t,x)I(t,x)] |p = -R| (t, 1). (6) 


Introduce the auxiliary function 


p(t) = a zO(x)L(x)V (t, x) I(t, x)dz 


and Lyapunov function F(t) = E(t) +ep(t). Here e € (0,1/K) is a parameter. 
On one hand, we have |p(t)| < KE(t) by using condition Hz and Cauchy 
inequality and F(t) is equivalent to E(t): 

(1— Ke)E(t) < F(t) < (14+ Ke)E(#). (7) 


The parameter 0 < € < 1/K ensures that F(t) is positive definite. On the 
other hand, we have 


< alt) 


= [L(1) + C()R?]J (t, 1)? — 22) — f xO" (x)|V (t, x)|? + 


zL (2)I@, 2) Par- lt), (8) 
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It follows from the above equality and assumption H3 that 


d L(1) + C(1)R? 


Pt) < 5 I(t, 1)? — Et). (9) 


By differentiating F(t) and using (6), (7) and (9), we obtain 


d d d 
GF = GEO +eF 0) 
<- [r E0 +CWR)] 4 yp _ FW 


2 14+ Ke 


Finally, choosing € to ensure R — (e(L(1) + C(1)R?)) /2 > 0 and applying the 
comparison principle (see Section 3.1 of [33]) and (7), we get 


E(t) < Me™ E(0) 


with M = (1 + Ke)/(1 — Ke) and w = ¢/(1 — Ke). Therefore we obtain the 
exponential stability of the system (3) and complete the proof of the theorem. 


Applying the auxiliary function p(t), one can also prove an observability 
inequality. 


Theorem 2 Under the assumptions H,-H3, for any V?(-), I?(-) € L?(0, 1], there 

exist time T > 2K and a constant C7 such that the solution to the system 
C(x)Vi(t, x) = —Iz(t, z), TE (0, 1), 

L(x)Ii(t, £) = —Vz(t, £), t > 0, 


Vit,0)=0, I(t,1) = 0, (10) 
V(0,x) = V? (a), I(0,x) = 1°(a). 
satisfies the observability inequality 
f |V (t, 1)|? dt > C-E(0). (11) 
0 
Proof: It follows from (6) that the energy of the system (10) satisfies 
d 1 
SEU) = V(t, 2) a)l = 0. (12) 


Similarly, from (8) and (10) we have 


2200) 
= C(1)|V(t, 1)? — 2E) =j E O E 
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which implies that 


with the help of assumption H3. Integrating the above inequality from 0 to 
7 and using (12) and |p(t)| < K E(t), we get the observability inequality (11) 
with the constant C, = 2(T — 2K)/C (1). Therefore we complete the proof of 
the theorem. 

For convenience in the convergence analysis, we rewrite the exponential 
stability of the system (3) in the language of semigroup of bounded linear 
operators. 


2 
Remark 1 Let H = |e (0, 1] be the state space and define the inner product on H 


by: (m, im) 4 = fy C(a)V(«)V(a) + L(x) I(a)I(a) de, Ym = (V, I), ñ = (V, Ù € H. 
Define the system operator A on H through: 


p(z) ol daa) 
a(o) = a: grle)” on 
D(a) = { (P03) emo DP: 9(0) =0, pa) = Raa). aay 


Thus, the system (3) can be transformed into abstract Cauchy problem 


d (V(t, x) V(t, x) V(0, £) V? (x) 

= =A , =| o cH. 

dt \ I(t, x) I(t, x) I(0, x) I(x) 
A routine method as given in Theorem 3.1.11 of [34] can be applied to show that the 
operator A generates a contractive semigroup T(t) on H. This means that the system 


(3) has an unique solution for any V°(-), 1°(-) € L?(0,1). Theorem 1 further implies 
that the semigroup T(t) is exponentially stable under the assumptions H,-Hs3, i.e. 


IT) la < Me", 


here ||- |g is the norm induced by the inner on H. 


3 Uniform exponential stability and uniform 
observability 


Firstly, we introduce the semi-discretization scheme for the system (3). For 
this purpose, let N € N be a positive integer and h = 1/(N +1) mesh size. 
Insert N + 2 points and N +1 points, denoted by y; = ih (i = 0,1,--- ,N+1) 
and x; = (j + 1/2)h (j = 0,1,---,N) respectively, in the domain [0,1]. If 
let f; be the value of any continuous function f(x) at the node y; = ih (i = 
0,1,--- ,N +1), then the notations 


fizi- fj fi 
h 


fit fi 


ôsfj = 2 


bi fj = 
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denote the central difference operator of f(x) and the average operator of 
f(x) at the node x;, respectively. Inspired by the works of [21]-[26], we propose 
the following semi-discretization scheme for (3) 


Cy G52 VE) = l(t), Ly G11; (t) = —deVj(), J = 0,1,---N, 
Vo(t) =0, Vivyi(t) = ca ), 
V(O, y) = V° lyi) = V; I0, y:i) ~ I (y i) = I, 4=0,1,---N+1. 

(15) 
We call this semi-discretization scheme as average central-difference method 
since the average operator and central-difference operator are applied for the 
temporal derivative and spatial derivative, respectively. 
The energy of discrete system (15) is 


2 
1 | + Lj 


2 
ò, 1 (8) | l 
which is discrete counterpart of the energy E(t). 


Definition 1 If there exist two constants M and w independent of t and h such that 
En(t) < Me™® Ep (0), 


then we call that the system (15) is uniform exponentially stable. 


The state space of the discrete system (15) is Hy = [RN+]? with the 
inner product 


N 
((pns dh); (Uns Yn) yw =h >> leo + Ljôqjðzv] , 
j=0 


for all (pn, qa), (un, Un) € Hn. Here pnr = (p1; ,Pn41), qh = (Gos*** GN) € 
RNI (un, un) and (pnr, qn) are the same type of vectors. Moreover, set artifi- 
cially k = R7}, pp = uo = 0, qN+1 = kpn+1, and vy41 = kuyn+1, to unify the 
notations of Ô Pj and 6,u,; etc. When R = 0, i.e. qu+1 = 0, the space Hy is 
still meaningful and denoted by H9, in this case. 

To verify the uniform exponential stability of the systems (15), we will 
follow every step of the proof of Theorem 1. However, we need the following 
lemma additionally. 


Lemma 1 Let {u hh, {vi 5" and {wah be sequences consisting of real 


numbers, then we have 


N 
2 (u41 — Us) (Vi41 + V4) (wi41 + Wi) 


ha 
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1 N 
+3 2 (uiţi — ui) (Vit — vi) (wi+ı — wi) 


N 
1 
+I 2 uiti + Ui) (Viti — Vi)(Wi41 + Wi) 


1 
+3 S (uipi + ui) (vipi + vi) (Wig — wi) 
i=0 


= UN4+1UN41WN41 — UOVOWO 


Proof: Extracting the factors from the first two terms and the remaining 
terms respectively and by simple algebra operations, we obtain 


N 


1 
Left = z X [lui = Ui) (Vig Wid + wivi) + (ui+1 + Us) (Vig Witi = wii). 
i=0 


Breaking the brackets in right hand side of the identity above and eliminating 
the cross terms, one has 


N 


1 
5 Se Uipi — Ui) (Vig Wii + Wii) + (uipi + ui) (vii Wit — wivi)] 
i=0 


2. 


= X (Ui41Vi41Wi41 a uiviwi) = UN+1UN41WN+41 — UDUOWO- 


We complete the proof of the lemma by the identities above. 


Theorem 3 Under the assumptions Hı and Hə, the semi-discretized system (15) is 
uniform exponentially stable. 


Proof: Firstly, differentiating the energy Ep (t) with respect to time t along 
the solution to (15), we have 


d 
gO 


N 
= -hY [55 V; (05215 (2) + ô 1)(O52Vj 0) 


j=0 
N 
= -5 in = WO) +50) 
j=0 
1 N 
$5 Un) - FOV) + 40) 
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j=0 
= —[Vw4i(t)In+1(t) — Vo(t)Lo(t)]. (16) 

This means that q 
ge = —R|In y (t). (17) 


Secondly, assume that £ is a parameter and define Lyapunov functions by 
Fi,(t) = En(t) + epn(t). Here p(t) are auxiliary functions given by pp(t) = 
—h am C5L561 4561 V;(t)d1;(t). On one hand, it is easy to see that |p, (t)| < 
KE,(t). This implies that Lyapunov functions L;,(t) are equivalent to the 
discrete energy Fp(t), i.e. 


(1—eK)En(t) < F,(t) < (1+ eK)E,(t). (18) 


The parameter 0 < € < 1/K ensures that Lyapunov functions F(t) are 
positive definite. 

On the other hand, differentiating the auxiliary function p;,(¢) and applying 
(15), we derive 


N N 
EmO = WS Ogu GOS) AAY Ljðgvðg VOV) 
l i j=0 j=0 
E XO Lilujaa + yat) GAl) + E) 
j=0 
N 
+ XO Oluja + ya) (Vjat) — VEVE) + VE). (19) 
j=0 


N 
TE Litua tyre) — HOE a) + HO) 


j=0 
N N 
h h I t)|? 
= 2 aO- HOP — 2 LDO + Ly MEPL (ao) 
j=0 j=0 
and 
1 N 
EE Goa +s) Viv) — Vi) Wp (8) + VO) 
j=0 
N N 
h h 2 a, Wan? 
=z 26 Vlt) — V(t)? — > 2 54Vj(0)| + Ora) 
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Combining (19) and (20)-(21), we get 


Sonlt) < alina OP- Ent), (22) 


in which a = 1/2 (Ly + CyR?). Note that the inequality (22) is a perfect 
counterpart of inequality (8). 

Finally, differentiating Lyapunov function F),(¢) and using (17)-(18) and 
(22), we have 


Choosing £ to ensure R — ea > 0 and applying the comparison principle and 
(18), we get 


En (t) < Me“ E, (0). 


M and w are the same as in the proof of Theorem 1. Therefore we obtain 
uniform exponential stability of the system (15) and complete the proof of the 
theorem. 

Similarly, using the auxiliary function p(t), we can also prove the discrete 
versions of observability inequality 11 hold uniformly with respect to step h. 
This is the main content of the following theorem. 


Theorem 4 Propose the semi-discrentized scheme of the system (10) as follows 


Vo(t) =0, In+i(t) = 0, (23) 


Then under the assumptions H,-H2, for any (VP, oo Vaart, 19, oo ry (a He 
there exist time T > 2K and a constant Cr independent of h such that the solution 
to the system (23) satisfies the uniform observability inequality 


OAOE (24) 


Proof: It follows from (23) and (17) that 


Í BO = -R Inn] = 0. (25) 


Similarly, from (23) and (19)-(21) we have 


d C(1) 


qe = Viv I? — E,(t). 
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Integrating the above inequality from 0 to 7 and using (25) and |pp(t)| < 
KE; (t), we get the uniform observability inequalities (24) with the constant 
C, = 2(r—-2K)/Cn 41. Therefore we complete the proof of the theorem. 

At the end of this section, we explain the mechanism of boundary feedback 
control of our numerical scheme how to uniformly exponentially stabilizing the 
discrete system (15). 


Remark 2 Let Vp(t) = (Vi(t),--- ,Vn4i(t))' and In(t) = (Io(t),--- ,In(t))' be 
the unknown variables of (15). We solve Vi (t) from (15) by letting j = 0 in the first 
equation of (15) since Vo(t) = 0. Using V1 (t), we can solve V3(t) by letting j = 1 in 
the first equation of (15). Repeating this process, we can separate all components of 
Vi (t). But we obtain I; (t) by the converse process in view of In 41(t) = kVy41 (t). 
That is the last component Iy (t) is the starting point of solving I; (t) from the second 
equation of (15). Thus we obtain the equivalent form of the discrete system (15) 


ORAON = (Bil, (t),---, Byailn(t), DoVa(t),---, DNVa(t))', (26) 
in which 
ByIn(t) = 2) 
del; (t) 


By+1In(t) = -2 


Bnailn(t) = TOF BnIn(t), 
ÔxVyn(t 
DnVn(t) = -2 n Fa kByyiln(t), 
N 
Dit) = s Bea), i=0,1, e, N1 


a 
satisfy the recursion relations and are all mappings from RHI toR. 

If you track the feedback control kVy +1(t), which corresponds to kBy+,Ip(t), in 
the dynamical system (26), you can find out that it firstly enters into the channel 
I(t) and then every channel of the system (26) one by one. This is caused by 
the average operator for the time derivative. Without average operator, the semi- 
discretization scheme (15) degenerates into the classical central difference scheme 
and the feedback control only appears one channel. Furthermore, it has pointed out 
that there is no uniform exponential stability for this approximating scheme in [21]. 
More information is given in section 6. This is the main mechanism of boundary 
feedback control uniformly exponentially stabilizes the discrete system. This has the 
same effect with the mized finite element method in [18]. 


4 Convergence analysis 


We will show that the solution to the system (15) converges to the corre- 
sponding solution of the system (2) in the sense of Trotter-Kato. Since the 
convergence analysis for (3.8) and (10) is the special case of R = 0, we skip it. 

For every n = 1,2,--- , there exist bounded linear operators Pp : X > Xn 
and En : Xn > X satisfying 
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(A,) There exist two positive constants M, and Mə such that ||E,,|| < Mı 
and ||P,|| < Mə, 

(Ag) |En Phe — x|| > 0 as n > œ for all x € X, 

(A3) PaEn = In, where I, is the identity operator on Xn. 

The notation B € G(M,w, X) with M > 1 and w € R, means that B is 
the infinitesimal generator of a Co-semigroup S(t), t > 0, satisfying || S) || < 
Met. The Trotter-Kato Theorem for approximating a linear Co-semigroup 
S(t) on a Banach space X is as follows. 


Theorem 5 (Trotter-Kato [35]). Assume that (A1) and (A3) are satisfied. Let 
B resp. Bn be in G(M,w,X) resp. in G(M,w, Xn) and let S(t) and Sn(t) be the 
semigroups generated by B and Bn on Banach spaces X and Xn, respectively. Then 
the following statements are equivalent 

(a) There exists a Ao E€ p(B) A NZL: p(Bn) such that, for all x € X, 


|En(Ao — Bn) Pn — (Ao — B)~*2|| > 0, as n> oo. (27) 
(b) For every x E€ X andt > 0, 
\|EnSn(t)Pn2 — S(t)z|| +0, as n— oœ (28) 


uniformly on bounded t-intervals. 


Note that the assumptions B, E€ G(M,w, Xn), or equivalently ||S,,(t)||n < 
Met, n = 1,2,---, usually is called the stability property of the approx- 
imations, whereas statement (a) is called the consistency property of the 
approximations. However, one may face some major difficulties when one 
applies Theorem 5 to perform convergence analysis. The most difficult one is 
how to verify the consistency property (a). The following property, which can 
replace (a) by a condition involving convergence of the operators B,, to B in 
some sense, is useful in this part [35]. 


Proposition 1 Assume that the assumptions of Theorem 5 are satisfied. Then 
statement (a) of Theorem 5 is equivalent to (A2) and the following two statements: 
(C1) There exists a subset D C D(B) such that D = X and (oI — B) tD = X 
for a ào >w. 
(C2) For allu € D there exists a sequence (Un)nen with Un E€ D(Bn) such that 
lim Enu=u, lim EnBnū = Bu. (29) 
n— o0 


n— oo 


Now we give the convergence analysis of our systems. In light (26), the 
approximating operators Ay are obviously defined by, for any (prgn) € Hn, 


m Bidn Dopn 
An Ge = ( A , with Bqn = ; , Dpr = = |. 60) 
Dpn i 
Byiid DnPn 
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By the same operations as in (16), one has 


N 
“i K N j=0 i a 


This means that Ayn € G(1,0, Hy). However it follows from of Remark 1 that 
A generates some contractive semigroup, i.e. A € G(1,0, H). This shows that 
the discrete scheme of (15) is stable. 

Let xs be the characteristic function of the set S and define the extension 
operators Ey : Hy > H: 


N 
En (>) = Dino (84 Pi)X yiyi] 
dh i= (52 G)X(yisyisal 
Choose the dense subset D = D(A) ()(C?(0, 1])? of H. For any (u(-),v(-))" € 


D, set U = (u(y1), oe u(yn4i)) | v= (v(yo), acta v(yn))' and U = 
(u,0)". By u(yo) = u(0) = =0 and v(1) = v(yn41) = ku(yn+1) = ku(1), it is 


easy to see 
N ulyi +u(yi 
EU E (7) = Di=o | Xlyiviti] 
NU=LN\_}] = N cf lesa bates 
. Lio —— | X[yiyi+ı] 
and 
N V(Yi —v(yi 
EyAnU = E Bu = Lixo oe Torta Xlyiyita] 
NANT EN | Ca) T | 55 (tun-u) 


i=0 hLi X[yiyita] 


N v(yi +v(yi N 
ys SE Xivani] — Limo UF) X bys svisal 


Ce 1) — ule) bulge) — e 
(yit1) — 


v(x) + u(y) = U2 ives dl 
a 
2 Licolv (E? 


41) Wiri — T) +u (EF 

1) (Yiz — 2) +V (E 
in which the mean value Theorem is applied and u(x), v(x) € C!1[0, 1] implies 
that |u’(&)|, Jw (Et), v (E)| and |v’(€?,,)| are uniformly bounded with 


N u(yi41) tuys N 
o avisi) tulu) Xlyiyitil T Xio U(X) X i yi1] 
1 
2 


ie Yi) Xiu vi] 
)(x — Yi) Xluvi] ' 
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Uniformly exponentially stable approximation of transmission line 15 
respect to i = 0,1,--- , N. Let T be their common upper bound and so we have 
T u(x)\ |? TR IES, Xy, : 
le» (7) — ( ) z ig? Alvia] +0, as N > oo. 
v v(x) N 2 Ži=0 X[yisyital N 
Similarly, for Ey Ay — A we obtain 
Bwaw (3) -4 (u) 
N v(yi —v(yi vu! (x 
_ eo k tal si ets X[yisyital 
_ N uļyi —ulyi ul (x 
ae w H (a) meh Xlysyi+t] 
N [o (nessa) vu (x 
_ Vi=o a = a Xlysyi+t] 
E N [u (na) _ u (£) 
ae, | Li E sn] X[yiyiti] 
N  C(«)[v' (ni.1)—v(@)|4[C(x)—Ci]v' (x) 
5 i=o = CiC(a) Xlyeyi+a] 
E N  L(a«)[u! (n 1) —ule)]+E(x)- Lilu’ (x) 
Vi=o = LiL(a) Xiu vit] 
By the same idea as above, we can show that Ey Ay 2 uz) converges 
U v(x) 


to zero as N + oo since C(x), L(x) € Ct[0,1] and u(x), v(x) € C?[0, 1]. 
Therefore, the statement (b) in Proposition 1 holds. 

Finally, construct the projecting operators Py : H > Hy in light of the 
expressions of the extensions Ey by 


Py (Wa) = (raa) 


with 


and 


( 
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Hiti 
u(t) = no | v(x)dz — Ij41u(t), 7 =0, 1, ---, N-1. 
Yj 
It is easy to show that Ey and Py satisfy (A1) and (A3). 


To prove that Ey and Py satisfy (A2), we firstly assume (u(x), v(x)) € 
D(A). With this assumption we have 


N i+ 
En Py e — p! Dizo y ! u(s)dEX ly: yi1] = Tico UF) Xtyswer 
y D ie v(x) d2X[y;,y:41] Jio V(OP )X iu yi 1] 


and 


u(z)) (a) _ [ERO -olx 5 9, ae NY soe 
wry (a) À EO a E ee 


here 6“ and 6? are chosen such that tie u(x)dx = u(0%)h and Ti v(x)dz = 
v(0?)h when the mean value theorem is applied, and the continuity of u(x) and 
v(x) on the internal [0, 1] is applied in the last step. Thus combing this result 
and the density of D(A) in the state space H, we obtain (A2). Moreover, let 
(V°, 1°) € H be the initial value of (10) and set (Vp (0), In(0))' = Py (V®, 1°) 
be the initial data of (15), then (A2) implies that (V;,(0), ,(0)) convergent to 
(V°, 7°) in the sense of 


En(V;,(0), 1,(0))' + (V°, 1°)", as N > oo. 


In a word, we have completed the verification of (A;)-(A3) and (C1)-(C2) 
and this means that the solutions to the discrete systems (15) strongly converge 
to the solution of (2), i.e., as N + ov, 


levtvwrs eS _ T(t) e? l +0, Y ee) eH. (31) 
5 Return to the original system (1) 


Now, all results for the systems (3) and (15) are going to be translated into 
the system (1) and it’s semi-discretized systems, respectively. Recall that the 
state space H and the system operator A corresponding to (3) have defined in 
Remark 1. Similarly, we introduce the state space Ho and the system operator 
Ao for the system (1). The space Ho is [L?(0,1)]? with the inner product 
given by: for any (pC), a())7, (u(:),0(-))" € [Z2(0, DP 


(ED) = [O , ae, 
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Define the system operator Ao on Ho through 
(«)\ [7i (ce 
Ux _ dx \ C(x) 
Ao ee = da u(x) , (32) 
dz \ L(x) 


0, u(1) = Roy}. (33) 


D(Ao) = { e 

Thus, the system (1) can be transformed into abstract Cauchy problem 
d /(Q(t,x)\ _ Q(t, x) Q(0,z)\ _ (Q°(z) 

(Geen) =42(Geen)> Geos) = Gras) Fo 


dt \ (t, x) 
Furthermore, the relations V = Q/C and I = ¢/L determine a mapping 


W:Ho > H given by 
20) e Ho 


; ) CIMES 


0 Lig) o(x) 


I(x) 
The operators A, Ao and W have following basic properties. 


(3) =*(82) = 


Proposition 2 The operators A, Ag and Y satisfies: 
(1) The operator Y is an isometric isomorphism from Ho to H. 


(2)Let WU! be the inverse operator of V, then the operators A and Ag are similar, 


ie. Ag =U TAW. 
Proof: (1) Obviously, ¥ is a linear operator from Ho to H and invertible. 
ae» and also a linear operator from H to 


The inverse of Y is Y7! = 0 Lle) 
Ho. This implies that the operator Y is an isomorphism. Because the identity 


e ORA -e 
w € D(Ao) 


C(x) 
holds, we know that W is isometric. 
(x) = -1 
x € D(A) = D(Ao), we have U ie 


( 


2 
Ho 


2 


A 


(2) For any y 


aon (05) =- a) =G) 


which gives Ag = UAW. 
Now, we can give exponential stability result of the original system (1) 
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Theorem 6 Let the assumptions H,-H3 hold and the semigroup on the space Ho 
generated by the operator Ag be To(t), then To(t) is exponentially stable. 


Proof: Recall that the semigroup T(t) given in Remark 1 satisfies 
ITa < Me™ for some positive constants M and w. However, it fol- 
lows from Ao = W~1AW and similar semigroup theory in [36] that To(t) = 
W-!T(t)WV. Therefore, from (1) of Proposition 2 we have 


Zo()llzo = II TOYA = ITO < Me™, 


which means that To(t) is exponentially stable. 

Certainly, we can obtain Theorem 6 by the same method in Theorem 1. 
However, this indirect method for obtaining exponential stability of the orig- 
inal system (1) is important in studying the uniform exponential stability of 
it’s discrete systems. From semi-discretization scheme (15) and the relations 
V;(t) = Q,;(t)/C; and I;(t) = ¢;(t)/L;, we obtain semi-discretization scheme 
for (1) 


Cj Tô, 20) = -ő (4 - ; 
5 (t) \ _ Qi(t) Tm 
Lj 5,63 ie js 52 ( Get j=0,1,---N, (34) 
(t) = ) = ROn Lyp bN+ (t), 


2, (0, yi) ¥ ply) =, i =0, N1. 
The energy of discrete system (34) is 


DRED 


Now we show that the system (34) is uniform exponentially stable in the 
sense of Definition of 1. The state space of the discrete system (34) is Hon =: 


[R+] ? with the inner product 


2 


Cj 


h N 
Eorl) = 5), 
j=0 


N oes a gs 
va ees Bi. i a. v; 
(Pas Gn), (Uns On) on = DD CER + Lið, 7-0; z] l 


in which (Ph, qh), (Ùn, Ùn) € Hon,Po = Uo = 0, Gnai = K'PN+1, ÙN+1 = 
k/tn41, and k’ = ROR} LN are used. Let Vy be matrixes of order 2N +2 


given by 
1 1 1 1 
Wy = di soit E 
N ine { >, en’ Lo’ ah 
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which is well-defined since the assumption Hj is right. Using these matrixes, 
we define the linear operators from Hon to Hy by 


Ph Dh Dh 
= ~ |, V( 2") € Hon. 
(>) A 2) A ON 


Thus we can get the results similar to those of Proposition 2 and Theorem 6, 
which are main results of this paper. 


Proposition 3 For any fixed positive integer N, the operator Vy is an isometric 
isomorphism from Hon to Hy. 


Proof: Obviously, for a fixed positive integer N, Wy is an isomor- 
phism from Hon to Hy since it is invertible. The inverse of Vy is Uy = 
diag {co,--- ,cn, Lo, +++ , Ly}. Because the identity 

2 = 
_ Ph 
-|@) 


ae (2) 


holds, we know that Vy is isometric. 


2 N 


+h Lj 


j=0 


2 2 


Hyn j=0 Hon 


Theorem 7 Define operator Aon by the formulas Aon = Wy Ant. Then the 


abstract Cauchy problem 
d (Qn(t)\ _ Qn(t) 
dt ee = Aon ey a 


determined by the operator Aon is equivalent to the discrete system (34) without 
initial data. Furthermore, the discrete system (34) is uniform exponentially stable 
with respect to energy Eo; (t). 


Proof: It is easy to see that the discrete system (15) is equivalent to 
d /Vi(t)\ _ Va (t) 
ae (nin) =A (n) i 
from (26) and (30). Let Qa = (Qi(t),---,Qnai(t)), alt) = 


(do(t),--+ , dn(t)) and (Qn (t), bn(t))' be state variables of the discrete system 
(34). Then we have the identity 
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Substituting the identity above into (36), we obtain 


sity (Gh) = An (GUO) 


which is equivalent to the discrete system (34). This means that the abstract 
Cauchy problem corresponding to the discrete system (34) without initial data 


m (ott) on (Gala) = earen (SiG) 


That is to say Aon = Vp An Wy. Let Ton (t) and Ty (t) be semigroups gen- 
erated by Aon and Ay respectively, we have To(t) = V~!T(t)W. Therefore, 
from Proposition 3 and Theorem 15 we know 


Zon (lon = [IF Tw) Vllzon = Tw lla < Me™, 


which means that Ton(t) is uniform exponentially stable, i.e. the the expo- 
nential decay rates of Ton(t) are independent of N. But if you restate this 
result in the language of energy, you can obtain that the discrete systems (34) 
are uniform exponentially stable with respect to energy Eo; (t). 

Finally, the solution to the discrete system (34) is convergent to the one of 
the continuous system (1) in the meaning of Trotter-Kato Theorem. In fact, 
constructing the extensions Hon from Hoy to Ho by Eon = U1 En Uy and 
the projecting operators Pon from Ho to Hon by Pon = Uy EnV, we have 
the following convergence result. 


Theorem 8 As N > co, we have 


EonTon(t)Pon ey - Tolt) a) 


P +0. ¥ e € Ho. (38) 


Proof: From the relation To(t) = Y~'T(t)W of the semigroups, Proposi- 
tion 2 and Proposition 3 we have 


UL 
[EonTon(t)Pon — To(t i (x) | 
Ho 


= |W" [EnTn (t) Px — TH)]Y ew (K Yh. 


= |[[EnTw(t)Py — T(t oye 


The identities above and (31) imply that (38) holds and the proof of the 
theorem is finished. 
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6 Numerical simulations 


In this section, we show the effectiveness of our numerical approximating 
schemes (15) or (36) and (37) through some numerical experiments. Recall 


that Uy = diag {2,--- : i et , L} and Aon = Uy AnUy are given 
in last section. Because the operator Ay and Aon are similar, we only give 
eigenvalue distributions of Ay to analyze uniformly exponential stability of 
(15) or (36). For this purpose, we should express the operator Ay as matrix. 


Let G), = diag{0,--- ,0,1}, 


11 1 


ʻi -11 
belong to CN+)xW+1) | Set 


1/B} 0 1 (-kG;, Mg 
y = ee a7 and Oy = 5 (TAY 0 


Then (36) is equivalent to 
a ie (t) Vi (t) 
= Val 39 
in which Ay = PK Vn On is used. It is easy to see that ®y is corresponding 


to the average operator of time derivative of (15). If one replace Bp by identity 
operator, then the classical finite difference scheme of (3) is easily restored 


from (39), i.e 

d /Vi(t)\ _ Va (t) 

dt aA ee GH ey) 
If C(x) = L(x) = 1, the numerical approximating scheme (15) degenerates 
to the numerical approximating scheme (2.5) of [25]. The authors pointed out 
that the numerical approximating scheme (2.5) in [25, Lemma 5.1] is essentially 
similar to the mixed finite element scheme of [12, 18]. However, this is not 
accidental, mainly due to the intervention of the average operator. 

Now we explain the significance of the discrete scheme (39) or (15). We 
plot two figures in Figure 1 and Figure 2, respectively. Figure 1 depicts the 
maximal real parts of Ay and VyQy for N = 40:5: 400. Figure 2 depicts 
the distributions of the eigenvalues of Ay and VyQy with N = 500. 

We see that the real parts of the eigenvalues of Y yQy approach to zero and 
those of Ay approach to a negative number from both figures. In both figures, 
we take k = R = 1, C(x) = In(1 + x) and L(x) = e”. Numerical simulation 
results show that the classical finite difference scheme of (40) is not uniformly 
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Fig. 1 Maximal real parts of eigenvalues of the semi-discrete schemes 


exponentially stable. This conclusion is consistent with earlier research results 
of [21]. However, Figure 1 manifests (15) is perhaps uniformly exponentially 
stable and this is in accordance with theory results of section 3. 


7 Concluding remarks 


Transmission line is a basic structure of circuit and play important role in 
physics and engineering. This paper is devoted to uniformly exponentially 
stable approximations for the transmission line with varying capacity and 
inductance. This means that we study it from the viewpoints of the numerical 
approximating and control theory. It is well-known that there are many dis- 
cretization methods to discretize the spatial variables. It is nontrivial to pick 
one which preserves exponential stability among so many semi-discretization 
methods. On the other hand, if the capacity parameter and the inductance 
parameter are constant, many existing results can be applied and there is no 
any challenge. To bypass the troubles brought by variable coefficients, suitable 
similar transforms are introduced. Uniform exponential stability and uniform 
observability are then smoothly obtained based on the method for the wave 
equations. These results have potential applications in uniform controllability, 
the approximation of the control problem and state reconstruction etc. They 
deserve to be investigated at length in the further research. 
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x without the average operator 
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Fig. 2 Maximal real parts of eigenvalues of the semi-discrete scheme with N = 500. 
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